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Abstract 

We demonstrate that the vortex state in the Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) phase 
may be very different depending on the field orientation relative to the crystalline axes. We 
calculate numerically the upper critical field near the tricritical point taking into account the 
modulation of the order parameter along the magnetic field as well as the higher Landau levels. 
For s-wave superconductors with the anisotropy described by an elliptical Fermi surface we propose 
a general scheme of the analysis of the angular dependence of upper critical field at all temperatures 
on the basis of the exact solution for the order parameter. Our results show that the transitions 
(with tilting magnetic field) between different types of mixed states may be a salient feature of the 
FFLO phase. Moreover we discuss the reasons for the first-order phase transition into the FFLO 
state in the case of CeCoIns compound. 
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I. INTRODUCTION 



Recent experimental studies of the superconducting state of CeCoIn 5 (see [1] and ref- 
erences cited therein) provided evidences in favor of the Fulde-Ferrell-Larkin-Ovchinnikov 
(FFLO) phase existence in the high magnetic field region of the superconducting phase 
diagram. Originally j3, S] the nonuniform FFLO state has been predicted to exist in su- 
perconductors when the magnetic field is acting on the electron spins only (the case of the 
paramagnetic effect). Usually it is an orbital effect which is the most important and this 
makes difficult the experimental observation of the FFLO phase. Moreover the supercon- 
ductor must be in the clean limit because the electron scattering is detrimental for the 
FFLO phase [4]. The orbital effect may be weakened in heavy fermion superconductors or 
in quasi-2D superconductors when magnetic field is applied parallel to the superconducting 
planes. That is why in addition to the heavy fermion superconductor CeCoIn 5 , quasi-one 
and quasi-two-dimensional organic superconductors are considered as good candidates for 
the FFLO phase realization (J El- Recently evidences of the FFLO state have been revealed 
in organic quasi-2D superconductors A-(BETS) 2 FeCl 4 [7] and ft-(BEDT-TTFS) 2 Cu(NCS) 2 
A- 

In the framework of an isotropic model with s-wave pairing the critical field for the 
FFLO phase in the presence of the orbital effect has been calculated by Gruenberg and 
Gunther [9] . They demonstrated that the FFLO state may exist if the ratio of pure orbital 
effect H°2 b (0) and pure paramagnetic limit H p (0) is larger than 1.28, i.e. the Maki parameter 
% = V%H°2 b (0) j H p (0) is larger 1.8. The pure paramagnetic limit at T = can be estimated 
as Hp(0) = Aq/a/2/xb, where A is the BCS gap at T = and \±b is the Bohr magneton flo| . 
In [9|] the exact solution for the order parameter was described by an FFLO modulation along 
magnetic field and the zero Landau level function for the coordinates in the perpendicular 
plane. Further analysis revealed that the higher Landau level solutions (LLS) become 
relevant for large values of Maki parameter c^m > 9 and the H c2 (T) curve may present 
regions described by different LLS. These results obtained for an isotropic model are readily 
generalized for the case where the electron spectrum anisotropy is described by an elliptic 
Fermi surface [l3|] . In such a case the Maki parameter becomes angular dependent and the 
transitions between different LLS may occur with a change of orientation of the magnetic 
field. 
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It happens that for an adequate description of the FFLO state in real compounds, the 
form of the Fermi surface as well as the type of the superconducting pairing play a very 
important role because they determine t 



cumstance has been demonstrated 



15. 



le direction of the FFLO modulation. This cir- 
281 ] in the framework of a general phenomenological 
approach based on the modified Ginzburg-Landau (MGL) functional fl4]. This approach is 
adequate near the tricritical point (TCP) in the field-temperature phase diagram. At the 
TCP the three transition lines meet: the lines separating the normal metal, the uniform 
superconducting state and the FFLO state. Near the TCP the wave vector of FFLO mod- 
ulation is small and this situation may be described by the MGL functional. For the case 

n 

when the deviation of the Fermi surface from the elliptical form is small the method [15|] 
permits to calculate the critical field corresponding to different LLS. 

(Section [III) Unfortunately the approach [l5|] is limited to weak deviations from the elliptic 
Fermi surface. In Section [XT] we develop a numerical method for the calculation of the upper 
critical field applicable to any cases, using a general form of the solution for the order 
parameter as a superposition of the different LLS. Note that the single LLS is an exact 
solution for the order parameter only for isotropic or quasi-isotropic (elliptic Fermi surface) 
cases. Otherwise the order parameter is described by an infinite serie of the LLS. However 
there is usually some dominating LL n and the amplitudes of other LL rapidly decrease 
with an increase of \n — n |. Our analysis qualitatively confirm the conclusions of [15] and 
reveal the transitions between the FFLO states with different dominating no. The obtained 
results demonstrate that the FFLO state, depending on superconductor parameters and/or 
magnetic field orientation, may take the form of the higher LLS with or without a modulation 
along the magnetic field. The transitions betwen these states result in a very rich dependence 
of the transition temperature on the magnetic field orientation. 

(Section HTT1) The approach of Section [III based on MGL is adequate near the TCP. On 
the other hand the case of anisotropic superconductors with elliptic Fermi surface may be 
treated exactly at all temperatures. In section [TTT] we use the scaling transformation J]^] 
to obtain the solutions for the higher LLS. As an illustration we consider quasi-lD and 
quasi-2D superconductors. Note that the higher LLS in quasi-2D superconductors with the 
in-plane orientation of the magnetic field were studied by Shimahara [if]]. 

(Section HVl) The experiments on CeCoIn 5 show that at low temperature the supercon- 
ducting transition becomes of the first order [Tjj . In Section [IV] we argue that this may 



be explained by the magnetism associated with cerium cites and its interaction with the 
superconducting subsystem. In the framework of the proposed model in the mixed state the 
cerium polarization should strongly increase in the cores of vortices. This mechanism could 
contribute to the anomalously large form factor of the vortex structure observed in CeCoIn 5 



at low temperature 



22]. 



II. FORMATION OF THE DIFFERENT FFLO STATES UNDER THE INFLU- 
ENCE OF THE ORBITAL EFFECT 

In this section we provide a general numerical approach for the calculation of the upper 
critical field of the FFLO state. Keeping in mind CeCoIn 5 we will consider the case of the 
tetragonal symmetry. Usual GL functional contains only the first derivatives of the order 
parameter and therefore it may be transformed by simple scaling transformation of the 
coordinates to the isotropic form (with the corresponding renormalization of the magnetic 
field). That is why in the framework of the GL approach we may easily obtain the exact 
solution of the H c2 problem for any anisotropic superconductor - the order parameter is 



decribed by the n=0 LL function j29|]. It is possible to decribe the FFLO state near the 
TCP point by MGL which takes into accound the higher derivatives of the order parameter. 
In contrast to the GL functional the MGL functional can not be reduced to the isotropic 
case by scaling the coordinates (this is possible only for s-wave superconductivity with 
an elliptic Fermi surface 13]). Therefore in the most general case, after performing the 
scaling transformation which makes the part with the first derivatives isotropic, we have the 
following modified Ginzburg-Landau functional: 



J 7 =V*\a + L)V (1) 

with a(H, T) = ao(T — T CU (H)) where T CU (H) is the transition temperature into the uniform 
superconducting state. The differential operator has the general expression 

l = -9 E n ? + 7 ( E n ?) 2 + ^ + |i n - K) + f n '} + iK n 2 .}) . ( 2 ) 

3 3 

where Ilj = — idj + 27rA 7 -/$o (with j = x ) y ) z) ) $ is the flux quantum, and the anti- 
commutator {Oi, O2} = O1O2 + O2O1. Here the z axis is perpendicular to the basal plane. 



iptic one and/or the pairing different from s-wave type. In contrast with previous 
151 ] they are not assumed to be small. 



To simplify the discussion, we assume that the effective mass is isotropic (taking account of 
its anisotropy is detailed in Appendix). For the appearance of the FFLO state g must be 
positive. The coefficients e Z) e X) e describe the deviation of the actual Fermi surface from 
the eL 
work 

The transition temperature is given by T C (H) = T CU (H) — A m i n /a where A m i n is the 
smallest eigenvalue of the operator L. Within the coordinate frame (x',y',z') with the z' 
axis pointing in the direction of the field, the eigenfunctions of L can be looked for in the form 
\P = exjp(iq z z f )(j)(x f ,y f ) where q z is the FFLO modulation vector along the field direction. 
In the abscence of anisotropic fourth-order terms, (j) can be found exactly. It is one of the 
Landau levels <p n defined in the (x' \y f ) plane. The eigenvalues are then 

X s °(q z , n) = 7 [((2n + 1)£ 2 + £ ~ Qof ~ it] (3) 
where n is a positive integer, the magnetic length 



«- s V^- <4) 

and the modulation vector 

qo = \[f. (5) 



In this case, A min = —7^ = ~9 2 /^l with a degeneracy of solutions (q z ,n) which is lifted 
when anisotropy is present. In the general case, we diagonalize L in the subspace of functions 
ip qzi n = ex P(^z^ / ) ( / :) n(^ / , y') (see details in Appendix) in order to find the smallest eigenvalue 
X(q z ) which is then minimized with respect to q z to get A min . 

A previous work [15] showed that due to the effect of anisotropy three different types of 
solution for the FFLO state can be realized: (a) the maximum modulation occurs along the 
magnetic field with the zero Landau level state, (b) both modulation and higher Landau 
level state, (c) the highest possible Landau level and no modulation along the field (or 
a modulation with a very small wave- vector) . Moreover due to the specific form of the 
Fermi surface a variation of magnetic field orientation may provoke transitions between 
states with different Landau levels. However a single- level approximation was used to get 
analytical results for these solutions. Due to this approximation the analytical results were 
valid for magnetic field higher than H ^> ^o-y -• I n the present work we show using 




FIG. 1. Angle dependence of the transition temperature and of the corresponding FFLO state for 
parameters e x = e = 0, with (a) e z — O.I7, g 507^ 2 , (b) e z — O.I7, 5 157^ 2 , and (c) 
e z — — O.57, g — 47^ 2 . The vertical bars in the bottom plots show the mean square deviation of 
the Landau levels composing the state from the average LL. 

numerical calculations that taking into consideration the full set of Landau levels, the results 
qualitatively remain true for arbitrary values of the magnetic field H. 

We calculate the transition temperature and the corresponding FFLO state when the 
magnetic field is applied in the xz plane. Typical results are displayed in Fig. [I] where 
we have set e x = e = 0. The form of the FFLO solution depends only on the parameter 
ratios e z j^ and ^h\J^ (see e.g. expression (IA6I) of the operator L in the basis of LL). As 
illustrated in Fig. [H(a), the FFLO state can appear with the maximum modulation vector 
q z = qo and the n = LL when the field is along the z axis. In contrast, for H along the x 
axis, there is no longitudinal modulation and the solution is composed by higher LL, which 
results in transverse modulations of the order parameter. When the field is rotated from 

5 




FIG. 2. Angle dependence of the transition temperature and of the corresponding FFLO state for 
parameters (a) e = e z = 0, e x = O.57, g 407^ 2 , (b) e x = e z = 0, e = — O.57, g 4Qj£~^ 2 , and 
(c) e = e x = — O.37, e z = O.27, g = 207^ 2 . The vertical bars in the bottom plots show the mean 
square deviation of the LL composing the state from the average LL. 

one axis to the other, the state is transformed with a continuous variation of the FFLO 
modulation and a smooth evolution of its expansion over the LL. However for smaller values 
of €hQo or e z /j the variation with the field orientation can be discontinuous with jumps of 
the FFLO modulation vector (see Figs.[T](b) and[T](c)). The sharp transitions are manifested 
by bumps and kinks in the angle dependence of the transition temperature. Fig. E9(a) shows 
that, for other anisotropy parameters, the jump can occur between states separated by more 
than one LL. As expected from the condition of single- level approximation H 3> $o^y^ ° r 
equivalently ChQo yjlj ^ 1? the number of LL that contribute significantly in the expansion 
of the FFLO state increases with the inverse of the field and/or the anisotropy (see Fig. [3]). 
The broadening of the expansion over the LL ends up in supressing the discontinuities. In 
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FIG. 3. Weight \c n \ 2 of the n-th Landau level (f n in the expansion of the FFLO state \I/ = 
ex.p (iq z z)J2 n c n ( fn at the field angle 9 = 7r/4, for parameters e x — e — 0, g — 507^ 2 , with (a) 
e z = -O.OI7, (b) e z = -O.I7, and (c) e z = -O.57. 



addition, as illustrated in Figs. [2](b) and E](c), the FFLO modulation and the transition 
temperature can vary non-monotonuously with the field angle. It is interesting to note that 
at the angles when the average LL is maximum the wave-vector of modulation is minimal 
(and vice versa) and it can even drop to zero (see Fig. 12(c)). At these regions the FFLO 
state corresponds to the highest LL states only. The experimental observation of such a 
non-trivial angular dependence of T c would be a strong evidence of the FFLO state. 
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III. HIGHER LANDAU LEVEL STATES IN THE FRAMEWORK OF THE MODEL 
OF EFFECTIVE MASS ANISOTROPY 



In this section we demonstrate how the higher LLS naturally appear in the exactly solvable 
model of the FFLO transition in a framework of anisotropic effective mass model. As it 
has been already noted in the case of the pure paramagnetic limit this model is reduced 
by a scaling transformation to the isotropic one with an arbitrary direction of the FFLO 
modulation [13]. In the presence of the orbital effect the situation is different and we consider 
here the uniaxial anisotropy (note that our results are readily generalized to the arbitrary 
anisotropy case). We are interested by a part of the Hamiltonian depending on the magnetic 
field T-Lorb +T-l>Pauii with 



1/3 ie „ A z 1 / d \ 2 1/3 ie 



^ orb 2m \dx ) 2m \dy ) 2m z \dz 



Upauii — I^bHg Z) (7) 

where we consider the effective mass m x = m y = m and the magnetic field H is in the 
xz plane making an angle 6 with the z axis. Our treatment can be readily generalized to 
the case of an anisotropic g factor [131 ]. The gauge of the vector potential A is chosen as 
A x = yH cos 0, A y = 0, A z = yH sin 9 and the spin quantization axis is along the magnetic 
field. 

Performing the scaling transformation z = z\f^ the orbital part becomes Jl^] 



1 / d ie TT \ 2 1 / d \ 2 1/9 ie TT I m . x 2 

ii H rn^ H 7 / H > / «it 



^ orb 2m \dx ) 2m \dy J 2m \dz f \j m z ' ) 5 ^ 

i.e. it corresponds to the isotropic metal with an effective mass m at the orbital magnetic 
field H = Hjcos 2 6 + ^- z sm 2 6 ( H z = H z and H x = H x \p~~~ )• The Pauli contribution 



may be written as 



T-t Pauii = ^bHcf z = VbH(7 z = /JjB h (T ^ (9) 

/cos 2 9 + ^ sin 2 9 



with the angular dependent effective Bohr magneton ji B (9) = /x^/^/cos 2 9 + ^- sin 2 9. 



In fact we have reduced the problem of the FFLO critical field calculation to that of the 
isotropic model with the field H and the effective Bohr magneton Jig (9) . The corresponding 
Maki parameter is c^m = V2H£ h (0)/H p (0) with, in our case, H^ b (0) that is determined by 
the effective mass m and then is the pure orbital field along the z axis, while H P (Q) = ^ . 
So the Maki parameter becomes angular dependent 

omW= W(0) = Q, 54 fM) . (10) 
A 0A /cos 2 # + ^sin 2 # 1 V dl Jt=t c , 



Remarkably in the later expression for cxm (0) enters only the slope of H c2 at the same angle 
9. As it was demonstrated in fll| for large values of the Maki parameter, % > 9, the critical 
FFLO field at low temperature is determined by higher LLS. In the case of a large quasi-2D 
anisotropy ^ 3> 1 this situation is realized when the Maki parameter is strongly increased 
for a field orientation near the xy plane. On the contrary for the quasi- ID anisotropy 
the Maki parameter is maximum for the field orientated along the z axis. 
The critical field may be numerically calculated from the formula 



ln (| 1 ) = ^ 27rRe E 



(11) 



where T co is the (zero field) critical temperature, uj n = 7rT(2n + l) are the Matsubara 
frequencies, L N are Laguerre polynomials, and 

T=T c0 



H (9) 12tt 2 V dT 



The LL number and the dimensionless vector of the FFLO modulation Q = hvpQ /2T co 
are chosen in a way to give the maximum critical field H (9). 

The calculated values of the upper critical field at T = K as a function of the critical field 
slope at T = T co are presented in Fig. BJ The LLS with n > appear at — ( ) > 

V / T=T c q 

18 (T/K). We see that with an increase of the slope the Landau level number n increases, 
while the FFLO modulation vector drops. For some slopes it occurs to be zero, and then 
the FFLO state is purely higher LLS. In the upper panel of Fig. [4] we observe the non- 
monotonous behavior of the upper critical field as a function of the slope (or the anisotropy 
2 ^ L ). With the increase of the slope the orbital effect is switched off and we approach the 
pure paramagnetic limit for the 3D case. However at T = K the transition into FFLO 




( dH 

FIG. 4. Upper panel: zero temperature critical field as a function of the initial slope I — it# 

\ ai j T=TcQ 

The transitions between the higher LLS are clearly seen. Lower panel: the FFLO modulation 
vector. 



state is a first order transition |12|] and then the calculated upper critical field should be the 
overcooling field of the normal phase. 

Large values of the Maki parameter suitable for the observation of these higher LLS 
states are mainly expected in strongly (quasi 2D or quasi ID) anisotropic systems. In 
such systems, the formation of the higher LLS may be clearly observed on the angular 
dependence of the critical field, which will reproduce the dependence on the initial slope 
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FIG. 5. The angular dependence of the upper critical field H C 2(0) at different temperature for the 
initial slope - (^ff^-)^ m = 17 (T/K). This case corresponds to the n — LL state. 



T=T c0 



displayed on Fig. [4] . In Fig. El such an angular dependence is presented for a maximum 

slope — ( dHc ^° ) ) = 17 (T/K), with a ratio of effective masses — = 100, below 

V ai J t=t c0 m 

the threshold of higher LLS formation. We see in Fig. [5] the standard behaviour inherent 
to the anisotropic mass model. The situation is very different in Fig. [6J where the slope 

— ( dHc ^° - ) = 60 (T/K) is well above the threshold. At low temperature the angular 

V /t=t c0 

dependence H c2 (9) clearly reveals the transition between the higher LLS, making the overall 
shape of the H c2 (9) curve very peculiar, and somewhat similar to the corresponding results 
of section II. 

Note that in isotropic systems the FFLO modulation vector Q is directed along the 
applied magnetic field. In anisotropic superconductor the FFLO modulation is described by 
exp (iQ (3^z f + - exp (iQ (^f^^^^ + V^TT^))' that is 

^ exp ( 7 ==JQ ( zcosO + — xsinO ) ) . Therefore the angle 9' that the direction of 

K V ^cos^+^sin 2 ^ V mz J J 

the FFLO modulation makes with the z axis is given by tan^ = — tan^. This means 

o J Vflz 
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FIG. 6. The angular dependence of the upper critical field H C 2{9) at different temperature for 



the initial slope 



dH c2 (90°) 
dT 



T=T c0 



60 (T/K) . At low temperature the transitions between the 



different LLS are responsible for the peculiar form of H C 2(0) dependence. 



that for quasi-2D anisotropy the FFLO modulation vector deviates from the field direction 
toward the z axis, while for the quasi- ID anisotropy it lies closer to the xy plane. 

The appearence of the higher LLS in quasi-2D superconductors when the magnetic field 
direction approach the xy plane is consistent with the prediction of such states in 2D super- 
conductors in tilted magnetic field 26]. 25]. 271]. 



IV. THE ORIGIN OF THE FIRST ORDER SUPERCONDUCTING TRANSI- 
TION IN CECOIN5 AND HIGH CONTRAST VORTEX PHASE 

In a magnetic field the superconducting transition in CeCoIn 5 becomes slightly first order 
below 0.3T CO for the field along the tetragonal z axis and below 0.4T CO for the field in the xy 
basal plane [17]. The change of the transition order occurs at a magnetic field lower than 
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that of the presumed FFLO transition. Another particularity of this compound is the field 
induced antiferromagnetic transition when the magnetic field is in the basal plane fisl . l^ . 
This antiferromagnetic phase exists only in the mixed state and basically in the same region, 
where the FFLO state is expected. Neutron scattering experiments [18| reveal a small value 
of the magnetic moment on cerium sites ~ 0.15/ig oriented along the tetragonal axis. We 
may expect that the normal state of CeCoIn 5 is very close to the magnetic instability of the 
itinerant type. Indeed the measurements of the magnetic susceptibility reveal its strong 
temperature increase at low temperature (at several K) f^ol ]. 

We propose to consider CeCoIn 5 as a system with two type of electrons, one mostly local- 
ized on cerium sites and responsible for the magnetism and the second strongly delocalized 
and responsible for superconducting properties. In fact a multiband picture for CeCoIn 5 
was already discussed by several authors - see for example [l| . 

In a magnetic field the electrons from the Ce band, which are polarized due to the 
exchange interaction, will create, in addition to Zeeman field, some internal field acting on 
the spins of the superconducting electrons. 

We may describe this situation by a simple Ginzburg-Landau functional introducing, in 
addition to the superconducting order parameter, the magnetic moment M of Ce sub-band: 

F(M, *) = -MH + A(T)M 2 + a(H + 1 M- H P (T)) |^| 2 + h - |^| 4 + 8F orh) (13) 

where H P (T) is the paramagnetic critical field and 5F or 5 describes the contribution of the 
orbital effect. The constant 7 describes the contribution of polarized Ce band electrons to 
the field acting on the spin of the superconducting electrons. Minimizing f fT3l) over M we 
find 

M = ' 1 14 

2A(T) V ; 

and finally substituting this expression into f fT3l) we obtain the pure superconducting func- 
tional 

* F *<*> - - um + a i H + "m] ~ H ' (T) ) '*' 2 + \( h ~ Sr> '*' 4 ■ (15) 

Here the role of the Ce band magnetization is the renormalization of the Zeeman field 
H H ^1 + 2A(T) ) an< ^ ^ e °^ coe ffi c i en t b of the |\^| 4 term b b — f^f)- Whatever 
sign of the exchange interaction 7, it decreases the coefficient of the |^| 4 term. With an 
increase of polarization at the normal state M/H = ^^y, it may even become negative. This 
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means that the superconducting transition becomes first order. We believe that namely such 
a situation is realized in CeCoIn 5 at low temperature, explaining the observed first order 
transition below (0.3 — 0.4) T co . 

Moreover the contribution to the magnetization from the Ce band (HID depends on the 
profile of the superconducting order parameter. In the vortex state the maxima of the mag- 
netization would be at the vortex core, where the superconducting order parameter vanishes. 
This circumstance may strongly increase the amplitude of the magnetic field modulation and 
then explain the anomalous behavior of the form factor of the vortex structure observed in 



CeCoIn 5 at low temperature 



21 



Previously Ichioka and Machida 



23| performed an extensive numerical analysis of the role 



of the Pauli paramagnetic effect in the context of quasi-classical Eilenberger theory. They 
demonstrated the modulation of electron spin susceptibility in the mixed state, increasing 
the electron magnetization in the vortex core and thus they have explained the anomalous 
behaviour of the form factor. In the recent paper [24], a similar problem has been treated 
by a variational method. Note that in these theories, only one electron band is implied, and 
therefore it is difficult to explain the occurrence of the first order superconducting transition. 
In some sense, our model provides an additional mechanism for the increase of the contrast 
of the magnetic field modulation. If we consider that the Ce band magnetism contributes to 
the temperature dependent susceptibility, while the band responsible for superconductivity 
gives the temperature independent contribution, then from the experimental data j^] we 
may roughly estimate that they are equally involved at low temperature (T ~ 2K). 



V. CONCLUSIONS 

In real compounds, the crystal structure plays a dominant role in determining the type 
of FFLO state. Of course it will also influence the vortex structure. The FFLO state may 
be characterised by an uni-dimensional modulation of the order parameter, and/or by the 
emergence of higher Landau level states. This is a crucial difference with superconductivity 
without FFLO state, where the crystal structure influence only on the type of the Abrikosov 
vortex lattice. The higher Landau level FFLO states should be realized in systems with 
strong uni-axial anisotropy and near in-plane orientation of the magnetic field. In such a 
case, the higher LL states should lead to an unusual angular dependence of H c2 . Finally 
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FIG. 7. Successive changes of coordinates: a scaling of the z-coordinate followed by a rotation 
around the y-axis. 



we propose a simple explanation of the first order transition at low temperature in CeCoIn 5 
based on the Ce-band magnetic contribution. 

Appendix A: Modified Ginzburg-Landau theory 

In the paramagnetic limit the MGL functional quadratic over \I/ is 



forth-order part of the IT, expansion. 
1. Rescaling of z coordinate 

The functional is invarient by the tranformations of the crystal symmetry group. In the 
tetragonal symmetry g x — g y — g ^ g z so, to recover an isotropic expression of the second- 
order part, one can rescale the z-coordinate as z = z \/ r m z /m x and the vector potential as 
A = (A 0x , A 0y , ^m x /m z A 0z ). Hence U Zo = ^m z /m x Il z so that £\ 9j^-% = 9 n j ( since 
gj oc TnJ 1 ). In doing so the magnetic field is transformed as 




(Al) 



3 



where a = a (T - T CU (H)), Uj 



idj + 27rA / /<I>o (with j = x,y,z) and L^Uj) is the 




(A2) 
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where H = V r x A and H = V ro x A . The angles that H and H have with the z-axis 
(see Fig. [Tj) are then related by the equality 



6 = arctan [ A / — tan 9 { 

m z 



rru 



(A3) 



Hence, after rescaling, the functional in the tetragonal symmetry has the general expression 



(A4) 



where the anti-commutator {0i,Oo} = O1O2 + O2O1. Note that in order to recover the 
functional used in previous work [151 ]. the term (27r/$ ) 2 [£(H% + Hy)+e x H%\ must be added 
in our expression. The latter term only shifts the energy by a constant so the solution for 
the order parameter is not modified. 



2. Expansion over Landau levels for the field applied in the xz plane 

In order to determine the transition temperature, one needs to find the eigenvalues of 
the operator L which is a polynomial of IF,. When the magnetic field is in the xz-plane, it 
is convenient to work in the rotated coordinate frame (x',y',z') where the z f axis points in 
the same direction as H (see Fig. [7j). By the change of coordinates the gradient operators 
are transformed as U x = cos 9H x r + sin 9U z f , U y = Ii y > , and Ii z = — sin 9H X / + cos 9U z f . 
Since the field H is along the z f axis, the operator Ii z > commutes with both Ii x > and U y f. 
So, with an adequate choice of gauge, the eigenfunctions of L can be looked for in the form 
^ = exp(iq z z f )(j)(x f ,y f ) where q z is the FFLO modulation vector along the field direction. 
In the abscence of the anisotropic forth-order terms, is a Landau level. Functions (p qzjn = 
exp(iq z z')(p n (x f ,?/), where p n are Landau levels, then form a natural basis over which to 
expand the solution in the anisotropic case. We use the orthonormal basis set composed by 
the states (p 2n = (?^) 2n (po/ \J (2n)\ and (p2 n +i = — i (ry1") 2n+l (po/y/ (2n + 1)!. Here is the 
operator of Landau-level creation defined as = ^§(IV ~~ ^ x >) where the magnetic length 

«» s S (A5 > 

and ip is the normalized lowest Landau level defined by r](po = 0. With v\ = ^75 (Fly + iH x ') : 
one can easily check for example that 777^ — 77^77 = 1 and Ft^i, + Fly = ^{2rfr] + 1). 
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After expressing the operator L as a function of rj and , the matrix elements L m ^ n = 
$ <p* qz , m L<p qx , n are found as 



_7_ 



"((2„ + 1 + A: 2 - A; 2 ) 2 - Ar 4 )^ n + L« B ] 



with 



= £hQ z and fc = 




(A6) 



(A7) 



They connect states which are separated by at most four levels. Within the above choice of 
basis set, the matrix is real symmetric 9-diagonal and the non-zero terms above the diagonal 
are given by the anisotropic contributions 

Ljtf = |l (4c 4 fc 4 + 12cV(2n + l)k 2 + 3s 4 (2n 2 + 2n + 1)) 

+ ^- (s 2 (3c 2 - 1) + 2s 2 (l + 3c 2 )(n 2 + n) + 2(1 + c 2 - 6cV)(2n + 1)A: 2 + 4cVfc 4 ) 
+ j- (c 2 (2n 2 + 2n - 1) + 2s 2 (2n + l)fc 2 ) , 



4 e) = (-l) n ^TTfcsin2^ 



.7 V^2 



s n 



+ 1) + V2c 2 A; 2 ) 



-(e) 



41 +3 = (-l)V(" + 3)!/n! 



V(n + 2)!/n!s 2 



ea;(w + 1) 
2^7 



" f ( (" + I) ^ + ^ ) + ( " C2(2n + 3) + (6C2 " 1)fc2 ) 



i. /l + 3cos2#. „. , 9 \ 

+ -^( (n + 1 - k cos 26> ] 

^2 7 V 2 / 



27 7 
fcsin2# 



V2 



7 ~ t" ~ 27. 



r(e) 



V(^ + 4)!/n!^ 



( )s c 

7 7 7 . 



(A8) 



with 5 = sin 9 and c cos 
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